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We present evidence of Kelvin excitations in space-time resolved spectra of numerical simulations 
of quantum turbulence. Kelvin waves are transverse and circularly polarized waves that propagate 
along quantized vortices, for which the restitutive force is the tension of the vortex line, and which 
play an important role iir theories of superfluid turbulence. We use the Gross-Pitaevskii equation 
to model quantum flows, letting an initial array of well-organized vortices develop iirto a turbulent 
bundle of intertwined vortex filaments. By achieving high spatial and temporal resolution we are 
able to calculate space-time resolved mass density and kinetic energy spectra. Evidence of Kelvin 
and sound waves is clear in both spectra. Identification of the waves allows us to extract the spatial 
spectrum of Kelvin waves, clarifying their role in the transfer of energy. 

turbulence in superfluids, whose main focus is to un¬ 
derstand nonlinear interactions and the resulting energy 
spectrum. Particular attention is given to the latter, as 
predictions for the energy spectrum differ. Using wave 
turbulence theory [2l| different teams of researchers ar¬ 
rived to different answers, one being that the energy spec¬ 
trum scales with wavenumber as ~ Q, and the 

other that it goes as ^ [T^ (it should also be 

noted that some authors claim that Kelvin waves play 
a neglible role in the spatial energy spectrum [1^). Re¬ 
cent numerical studies yield results compatible with a 
~ spectrum [23l - [^ . However, studies of wave dy¬ 

namics usually have a single or a few vortices, either us¬ 
ing vortex line dynamics simulations [ 2 ^, or solving the 
Gross-Pitaevskii equation (GPE) In the laboratory 
Kelvin waves have been detected in superfluid helium 
using submicron ice particles as tracers [20l| ; the analysis 
focused on selected reconnection events and the subse¬ 
quent emission of Kelvin waves. Finally, Kelvin waves 
and other vortex wave modes have also been studied and 
identified in BECs [2§ - l3Tj |. 

Our aim is to study Kelvin waves in quantum turbu¬ 
lence. As mentioned above, numerical studies of quan¬ 
tum turbulence have two main approaches. One is to sim¬ 
ulate quantum vortex lines dynamics [ 2 ^, ; in this ap¬ 

proach the velocity field outside the vortex lines is given 
by the Biot-Savart law and reconnection events are per¬ 
formed ad-hoc. The approach we use solves the GPE, 
an equation for the evolution of the wavefunction 'ip for 
a system of bosons. A Kolmogorov spectrum ^ 
was obtained in simulations of the GPE in [ij, Is^ - lssjl . 
We consider a flow that shares similarities with the von 
Karman flow generated in recent experiments with super¬ 
fluid helium [36| . We extract and save four-dimensional 
information to compute the spatio-temporal spectrum of 
quantum turbulence, where the presence of sound and 
Kelvin waves is evident. Previous studies using a simi¬ 
lar technique observed sound waves [s^, , and Kelvin 

waves in a single and straight vortex filament [ 2 ^. Our 


I. INTRODUCTION 

Quantum turbulence is the chaotic and erratic spatio- 
temporal behavior observed in superfluids and Bose- 
Einstein condensates (BECs) Qi- Its motion is char¬ 
acterized by the interaction between vortex filaments, 
where all vorticity is concentrated and which have quan¬ 
tized circulation Q. As superfluids have no viscosity, 
quantum turbulence has garnered much attention from 
the classical turbulence community, for it could provide 
insight into extremely developed turbulence. But it has 
also been a matter of debate how similar quantum and 
classical turbulence actually are. One of the charac- 
teristical features of turbulence, the existence of a Kol¬ 
mogorov energy spectrum, has been confirmed in super¬ 
fluids However, the physical mechanisms behind 

this spectrum are not completely understood. Moreover, 
other features, such as velocity statistics, appear to be 
different between classical and quantum turbulence iQ- 

An interesting property of classical and quantum vor¬ 
tex lines, first predicted by Lord Kelvin Q, is that waves 
can propagate along them when the line is subjected to 
helical deformations. These waves are known as Kelvin 
waves. In quantum fluids they are believed to play a 
crucial role in the turbulent energy cascade im , where 
energy is transfered from large to small scales and it is ul¬ 
timately dissipated by phonon emission [mil: while at 
large scales vortex interaction and reconnection flB | dom¬ 
inate the transfer of energy, at small scales Kelvin waves 
are believed to interact nonlinearly exciting fluctuations 
at even smaller scales. They pla y a role in various prob¬ 
lems in classical fluid dynamics [ 13 , but not in clas¬ 
sical turbulence, being the turnover time or the period of 
waves particular to that system (such as inertial waves 
in a rotating flow, or Alfven waves in a magnetofluid) 
the only relevant timescales [3, [l^. As a result, detec¬ 
tion of Kelvin waves in a disorganized superfluid flow is 
considered a signature of quantum turbulence SUS- 

There are several theoretical studies of Kelvin wave 
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study is done in a three-dimensional highly turbulent en¬ 
vironment, with a large number of vortices. 


II. THE GROSS-PITAEVSKII EQUATION 


The GPE describes the evolution of the wavefunction 
■i/i of a field of weakly interacting bosons of mass to, 

where g is proportional to the scattering length. By 
means of the Madelung transformation 





im4>{r,t) / h 

^ 5 


( 2 ) 


where p is the density of particles and the phase (j) can 
be associated with a velocity by v = Vd), one obtains a 
hydrodynamic description of the system [39j, 


dt 


+ V • (pv) = 0, 


9v g -n 



( 3 ) 

( 4 ) 


These equations are similar to the Euler equations for a 
classical and compressible barotropic fluid, except for the 
extra second term on the r.h.s. of Eq. ([H which is referred 
to as the “quantum pressure”. By solving the GPE we 
get the full three dimensional velocity and density fields, 
and compared to other methods [i^, we do not need to 
reconstruct them from the vortex tangle configuration. 

The hydrodynamic description allows us to define a 
quantity akin to the classical kinetic energy of a fluid, 
namely Ek = pv^ 12. Note that this is only a fraction of 
the total energy density, given by /2m +g\ijj\‘^/2. 

For simplicity, we will refer to the classical kinetic energy 
just as the kinetic energy. It can be further decomposed 
into an incompressible component and a compressible 
component E^, by decomposing the velocity field into ir- 
rotational and solenoidal components. The two compo¬ 
nents will be useful to discriminate between sound waves 
and other excitations in the fluid. More details on the 
energy decompositions are given below (see also [sMIdljl'). 

This system can have sound waves which follow the 
Bogoliubov dispersion relation [s^, 


ws(fc) = k\j c2 -I- (5) 

where the speed of sound is c = g\'il)\‘^/m and the coher¬ 

ence l eng th is ^ ^W'/{2m\ili\'^g) [^, , and Kelvin 
waves [3| which follow the dispersion relation 
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where a is the vortex core radius, and Kq and Ki are 
modified Bessel functions. The dispersion relation is 
quadratic in the small k limit, while for large k it is linear. 

To solve numerically Eq. ([T|) we use GHOST [i^, a 
highly parallel code which uses a pseudospectral method 
to compute spatial derivatives, fourth-order Runge- 
Kutta to compute time derivatives, and can solve PDFs 
in Gartesian periodic grids. We use 512^ grid points and 
the “2/3 rule” for de-aliasing. The speed of sound is 
c = 2 and the coherence length is ^ = 0.1/(8-\/2) in di¬ 
mensionless units in a three-dimensional box of length 
L = 2'k. Quantities are made dimensionless using char¬ 
acteristic length Lq, velocity Uq, and mean density po 
(see M)- These parameters result in an intervortex dis¬ 
tance £ such that ki = 2 'k/£ Ki 10 (see [H, H^), and in a 
vortex core radius a Ri 2^ as measured directly from the 
full width at half maximum of the mass density profile. 
Also, the quantum of circulation h/m in dimensionless 
units is given by 47rc^/-\/2 = O.OStt. 

As initial condition we use the Taylor-Green flow, 
which results in a set of vortex loops in two counter¬ 
rotating large scale eddies with turnover time of order 
unity, and whose geometry mimics the von Karman flow 
. The von Karman flow is used in recent experiments 
with two counter-rotating propellers such as SHREK (th e 
Superfluid High REynolds von Karman experiment) |36l| . 
and has been used in the past to measure Kolmogorov 
spectra in superfluids Q. For more information about 
the generation of the initial conditions see Appendix 
We let the simulation run for Ki 20 large-scale turnover 
times, so as to get good statistics on the slowest waves in 
the system. To also resolve the fastest waves, we use a 
very high output cadence, storing one output of the wave- 
function every half period of the fastest waves (i.e., the 
sound waves for the maximum wavenumber in the spatial 
domain). This very high temporal resolution allows us to 
properly calculate space and time resolved spectra. 


III. CHARACTERIZATION OF THE FLOW 

Figure [T] shows a three-dimensional rendering of the 
density field p. The visualization is done using the soft¬ 
ware VAPOR [i^, and only regions with low density 
(indicating topological defects associated with quantized 
vortices) are shown. The field evolves from a well ordered 
structure to a bundle of intertwined and structurally rich 
vortices. This evolution is associated with the develop¬ 
ment of reconnection events, which become prominent at 
t Ki 5 and start to subside after t Ki 10. The analysis 
below is done during this period, as afterwards too much 
of the total energy has decayed into phonons. Track¬ 
ing reconnection events and the subsequent generation of 
Kelvin waves in this flow is not an easy task. However, 
the identification of helical waves from spatial observa¬ 
tions is possible if single vortices or sing le reconnection 
events are followed in time (see, e.g., 

In Fig. [5] we present the evolution of the different com- 
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ponents of the total energy in the simulation. As men¬ 
tioned above, the classical kinetic energy density is de¬ 
fined as 


Ek = (7) 

The remaining components of the total energy density in 
the system are the so-called quantum energy 


2m2 


(vVp)^ 


and the potential (or internal) energy 


“ 2m2 




( 8 ) 

(9) 


A detailed analysis and derivation of each component of 
the total energy density can be found, e.g., in [4d| . 

From Fig. [5] different regimes can be identified in the 
evolution. In the first stage, up to t ~ 5, the incom¬ 
pressible kinetic energy oscillates around a mean value. 
Afterwards it decays and the other components of the 
energy grow, as reconnection of vortex lines takes place 
and the flow becomes more complex (see Fig. 1). After 
t ft! 10 the number of reconnection events in the flow sub¬ 
sides and the growth of the compressible, potential, and 
quantum energy becomes slower. During this process the 
total energy is conserved up to the sixth significant digit. 
Mass and momentum are also conserved. Snapshots of 
the velocity field and of the density between t ft 5 and 
ft 10, when the turbulent spectrum is more developed, 
are thus used for the spatio-temporal analysis. 


IV. LOOKING FOR NEEDLES IN A HAYSTACK 

Instantaneous flow visualization is insufficient to iden¬ 
tify and extract all the waves in a turbulent flow. In 
particular, in order to quantify their relevance in the en¬ 
ergy cascade, it is necessary to quantify their amplitudes 
as a function of frequency and wavenumber, i.e., to cal¬ 
culate space-time resolved spectra. As Kelvin waves are 
oscillations of the lines with p = 0, we should be able 
to identify their imprint in the mass spectrum, and we 
therefore study it first. Figure [2] shows the mass spec¬ 
trum p(fc,w), along with a close-up for low wavenum¬ 
bers, as well as a cut of p{k,uj) at k = 32. The dashed 
line in Fig. [2] corresponds to the Bogoliubov linear dis¬ 
persion relation of sound waves [Eq. ([5])] , while the solid 
line corresponds to that of Kelvin waves [Eq. ([HI)] . Sub¬ 
stantial power is concentrated along these curves. At low 
wavenumbers Kelvin waves are dominant, whereas sound 
waves become prominent as Kelvin waves begin to fade 
for /c > 50. At very high wavenumbers another accumu¬ 
lation of energy can be seen for ui ft 350, which is prob¬ 
ably due to quantum pressure effects inside the vortex 
core. Indeed, the spectrum of the quantum pressure also 
shows a bump at these wavenumbers (not shown), and 



FIG. 1. (Color online) Three-dimensional rendering of the 
mass density. The lines, corresponding to regions of low den¬ 
sity, are associated with quantized vortices. Typical length 
scales are indicated by the black bars (see text for descrip¬ 
tion). On the left we show a zoom into a single vortex; helical 
perturbations propagating along it can be identified. Due to 
the highly turbulent nature of the flow, identifying waves here 
is like looking for needles in a haystack. 


the frequencies are compatible with those pred icted for 
axisymmetric oscillations of the vortex core [3 [^. The 
presence of Kelvin waves at low wavenumbers is the most 
striking feature of the flow; these modes are probably ex¬ 
cited by deformation of vortex lines by the large-scale 
flow, and are hard to identify by simple inspection of the 
evolution of individual vortex lines. In Appendix |B1 we 
present a benchmark sutdy of the spatio-temporal spec¬ 
trum of mass in a flow where only linear Kelvin waves 
were excited, to verify the modes indentified in Fig. [5] 
correspond to these waves. 

To independently verify the presence of sound and 
Kelvin waves, and to separate the multiple branches of 
the dispersion relation, we now consider the spectrum 
of the compressible kinetic energy E^{k,uj), and of the 
incompressible kinetic energy E1(k,oj) (see Fig. jS]). A 
strong accumulation of energy around modes satisfy¬ 
ing the relation w = uisik) is evident in the compress¬ 
ible spectrum (Fig. |3J4). In the incompressible spectrum 
(Fig. 133) these excitations are negligible, but two new 
features are found. First, strong excitations can be ob¬ 
served for all wavenumbers at low frequencies. These 
excitations are compatible with sweeping of the vortex 
cores by the large-scale flow, i.e., the advection of small 
scale structures by the flow with a slow timescale asso¬ 
ciated with the turnover time. This is an important ef¬ 
fect in classical turbulence where it is responsible for the 
temporal decorrelation of modes [l^ [l^. It results in 
the smearing of the energy for all frequencies w = Urmsk 
and smaller (with r.m.s. velocity Urms ~ 0.5). Super¬ 
posed to these excitations, modes compatible with Kelvin 
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FIG. 2. (Color online) Space-time resolved mass spectrum, 
p(k,u)) normalized by p{k). A: The whole spectrum; the 
dashed (blue) line indicates the sound wave dispersion rela¬ 
tion tJs(fc), the solid (green) line indicates the Kelvin wave 
dispersion relation uJK{k). B: Close-up for small wavenum¬ 
bers; note the accumulation of power near the Kelvin wave 
modes. The range of wavenumbers with strong Kelvin exci¬ 
tations are marked as “inertial range”. C: A cut of p{k,u)) 
(also normalized) for k = 32, the marker indicates the Kelvin 
wave frequency. At low wavenumbers Kelvin waves are dom¬ 
inant, whereas sound waves become prominent after k ~ 50. 
At very high wavenumbers the quantum pressure inside the 
vortex core also leaves a trace in the spectrum. 


waves can still be observed for small wavenumbers. To 
further verify this, we calculated the decorrelation time 
Tfi of individual Fourier modes and compared it to the 
sweeping time Tg ~ l/(Urmsk) and to the wave period 
^ l/ujK{k) (Fig. [21). The decorrelation time fluctuates 
between Tg and for low wavenumbers and converges to¬ 
wards Tg for large wavenumbers; no other relevant time 
scales can be observed. Other collective vortex motions 
(e.g., Tkachenko waves [46l-l48l|) may also be present but 
masked by these two dominant timescales. A weak trace 
of these effects is also present in E^(k,uj). Kelvin waves 
do not contribute significant energy to the spectrum of 
classical turbulence, making the presence of energy in 
these modes a signature of quantum turbulence. 

The role of Kelvin waves in the dynamics of quantum 



FIG. 3. (Color online) A: Space-time resolved compressible 
kinetic energy spectrum, E^(k,uj), normalized by E%(k). B: 
Same for the incompressible kinetic energy spectrum El(k, uj). 
The dashed (blue) line corresponds to ujB(k), and the solid 
(green) line corresponds to uJK(k). The dash-dotted (black) 
line corresponds to sweeping, which excites all modes with 
frequency equal or smaller than u = Urmsk. El.(k,ui) is dom¬ 
inated by excitations around the sound wave dispersion rela¬ 
tion while El(k,uj) shows sweeping and Kelvin wave excita¬ 
tions. Inset: decorrelation time Td(k) of individual Fourier 
modes, compared with the sweeping time Ts (dash-dotted 
black line) and the Kelvin wave period Tu: (solid red line). 


turbulence is controversial. While at large scales the in¬ 
teraction of quantized vortices with the flow results in 
advection and reconnection of vortex lines, at scales com¬ 
parable to the intervortex distance £ Kelvin waves are 
believed to interact nonlinearly, exciting smaller fluctua¬ 
tions that eventually loose their energy to phonons. This 
transfer of energy by nonlinear interaction of wave modes 
can be described by wave turbulence theories, but cur¬ 
rent predictions differ on the shape of the energy spec¬ 
trum i®. The quantification of the amplitude of all 
Kelvin waves modes shown in Figs. ID and Fig. [3] can be 
used to shed some light on this problem. We thus ex¬ 
tract the modes centered around the dispersion relation 
given by Eq. © with a width of 2cr from the wave fre¬ 
quency, where the dispersion cr is estimated by fitting the 
peaks in p(ki,u!) such as the one shown in Fig. |2ld with 
a Gaussian. This subset of modes, which we call ^K^ 
can be used to compute the spectrum of the incompress¬ 
ible kinetic energy associated with Kelvin wave modes, 
Iqk HI- verified that the shape of 
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FIG. 4. {Color online) Spectrum of the incompressible kinetic 
energy associated with Kelvin wave modes, El{k,uj) dcj, 
where fix are the modes neighboring u)K{k). The two predic¬ 
tions for the Kelvin wave spectrum are shown as references. 
The shaded area corresponds to the region with strong Kelvin 
excitations identihed as “inertial range” in Fig. [2]B 


the spectrum is not very sensitive to the width in units 
of a used to define the subset fix, as long as it is not 
too large so as to include, e.g., modes associated with 
sound waves. In this spectrum, the candidate for an in¬ 
ertial range is observed for fc > fcf, and its width is in 
agreement with the wavenumbers for which we observe 
Kelvin-like excitations in Fig. As the scale separation 
is limited, we do not attempt to fit the data in Fig. |4] and 
only show the two theoretical predictions for the scaling 
as references. 


V. CONCLUSIONS 

Kelvin wave turbulence is an inherently quantum 
regime of superfluids. Proper quantification of Kelvin 
waves is important to understand the differences between 
classical and quantum turbulence. We presented direct 
evidence of the presence of Kelvin waves in numerical 
simulations of quantum turbulence using the GPE. By 
looking at the space-time resolved mass density spec¬ 
trum, we showed that Kelvin waves play the dominant 
role at scales comparable to the intervortex distance, 
while sound waves are excited at smaller scales. Further¬ 
more, the kinetic energy spectrum confirms the presence 
of both waves in the flow. The presence of these waves 
in the spectrum can be considered a quantum signature 
that distinguishes this flow from classical turbulence. 

Modes satisfying the wave dispersion relations excite a 
continuous region of the spectra, with broadening around 
the theoretical dispersion curves. This suggest that waves 
interact nonlinearly, as broadening of the dispersion re¬ 
lation in turbulent flows is often the result of nonlinear 
coupling and energy transfer. However, the modes ex¬ 
cited are still close to the linear dispersion relations, an 
unexpected result as excitations in the turbulent regime 
do not necessarily have small amplitudes. Finally, extrac¬ 
tion of the spectrum of the incompressible kinetic energy 
of modes along the dispersion relation of Kelvin waves 
results in a spectrum that is compatible with predictions 


of weak turbulence theories for superfluids. 


Appendix A: Preparation of the initial conditions 


The Taylor-Green initial conditions are generated by 
preparing a wavefunction ip whose associated velocity 
field is a Taylor-Green flow [s^, given by 

Vx{x, y, z) = sin(a;) cos{y) cos(z), 

Vy{x, y,z) = — cos{x) sm{y) cos{z), 

Vz{x,y,z) = 0. 

The Taylor-Green flow contained in a periodic box has 
properties that mimic the von Karman flow driven by two 
counter-rotating impellers, and has been extensively used 
in simulations to compare with the experimental flow. 
The von Karman flow is commonly used in the labora¬ 
tory to study turbulence, including dynamo experiments 
with conducting flows, experiments to study Lagrangian 
particles, and superfluid turbulence [s^. 

The process to generate the initial wavefunction is de¬ 
scribed in great detail in [s^. Here we give a brief, but 
nonetheless complete, presentation. The Taylor-Green 
flow can be described by the Clebsch potentials 

A(x, y, z) = cos(x)v'2| cos(z)|, 
lJ-{x,y,z) = cos{y)y'2\ cos( 2 ;)| sgn[cos(z)], 

which verify V x v = VA x Vy. Our use for them is 
that they can map a point in the (A, /i) plane to a line 
in three-dimensional real space. Now, instead of hav¬ 
ing to assemble directly a three-dimensional wavefunc¬ 
tion whose nodal lines match vortex lines of the velocity 
field V, we can pick instead a complex field (^(A, fi) which 
has a zero (a defect) at the point {Xd,yd)- Then, the 
three-dimensional wavefunction 


'fix, y, z) = (piXix, y, z), /r(x, y, z)), (Al) 


will be equal to zero along the line(s) defined by 
X{x,y,z) = Xd and y{x,y,z) = yd- This ensures that 
the defects of ip match vortex lines of v, as desired. 

Our choice for (p is 


(/)(A,y) =(j)e ( A - 


X (pe 








with 


</'e(A,/r) 


(A -I- iy.) 


tanh 




As (pe has one simple zero at the origin, (p will have four 
simple zeros in the region [0,7r] x [0,7r], resulting in four 
nodal lines in three-dimensions. But as the circulation of 
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FIG. 5. Time evolution of the incompressible kinetic energy 
E%, of the compressible kinetic energy E%, of the potential 
energy Ep, and of the quantum energy Eg. 


Appendix B: Analysis for helical perturbations 

As a reference, and for comparison with the spectra 
presented in Sec. IIVI for the Taylor-Green initial condi¬ 
tions, we present here the space-time resolved spectrum 
for a system of only four straight vortices at rest, per¬ 
turbed with small helical perturbations so as to excite 
linear Kelvin waves. The spatial resolution and parame¬ 
ters of the simulation were the same as in the simulation 
in the main text. Initial conditions were prepared as in 
[ 2 ^, and the vortices were perturbed with a superposi¬ 
tion of small helical displacements between wavenumbers 
k = 3 and fc = 30 (see also (H). The resulting space- 
time resolved mass density spectrum is shown in Fig. 
The perturbation excites strong sound waves, but modes 
compatible with the dispersion relation of Kelvin waves 
can be identified in the range of wavenumbers excited by 
the small initial perturbation. Moreover, the dispersion 
relation of these modes is compatible with the one ob¬ 
served in the main text for the fully developed turbulent 
flow. 


the velocity v must match the circulation generated by 
the vortex lines (whose circulation is in turn quantized 
and fixed by the parameters of the simulation), we must 
change the multiplicity of the nodal lines to match both 
values. Using that the total circulation of v is F = 8, 
that the quantum of circulation is dira (with a = h/2m), 
and that we want multiples of four nodal lines, we write 

(A2) 



where the brackets {} denote the integer part. 

This wavefunction is then evolved in time for over 30 
turnover times under the dynamics of the Advective Real 
Ginzburg-Landau Equation (ARGLE), 


dip 

'm 





—zv • Vip 



(A3) 


FIG. 6. Space-time resolved mass density spectrum p{k, uj) for 
a simulation with four straight vortices perturbed with small 
helical displacements between fc = 3 and fc = 30. The dashed 
(blue) line indicates the sound wave dispersion relation a;s(fc), 
and the solid (green) line indicates the Kelvin wave dispersion 
relation ujK{k). 
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